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Abstract. Using hydrodynamic computer codes, we 
study the possible patterns of relativistic expansion of 
an enormous pair-electromagnetic-pulse (P.E.M. pulse); a 
hot, high density plasma composed of photons, electron- 
positron pairs and baryons deposited near a charged black 
hole (EMBH). On the bases of baryon-loading and energy 
conservation, we study the bulk Lorentz factor of expan- 
sion of the P.E.M. pulse by both numerical and analytical 
methods. 
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In the paper by Preparata, Rufhni & Xue (1998), the 
"dyadosphere" is defined as the region outside the horizon 
of a EMBH where the electric field exceeds the critical 
value for e"'"e~ pair production. In Reissner-Nordstrom 
EMBHs, the horizon radius is expressed as 
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The outer limit of the dyadosphere is defined as the radius 
Tds at which the electric field of the EMBH equals this 
critical field 
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The total energy of pairs, converted from the static electric 
energy, deposited within a dyadosphere is then 



2 r+ rds \rds 
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In Wilson (1975, 1977) a black hole charge of the order 
10% was formed. Thus, we henceforth assume a black hole 
charge Q — O.lQmax, Qmax = VGM for our detailed nu- 
merical calculations. The range of energy is of interest as 
a possible gamma-ray burst source. 

In order to model the radially resolved evolution of the 
energy deposited within the e"'"e~-pair and photon plasma 
fluid created in the dyadosphere of EMBH, we need to 



discuss the relativistic hydrodynamic equations describing 
such evolution. 

The metric for a Reissner-Nordstrom black hole is 



~gtt{r)dH + grr{r)d^r + r^d^O + sin^ ed^<f) , (4) 
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wYieie gtt{r) = -grr{r) = - 1- 

We assume the plasma fluid of e+e^-pairs, pho- 
tons and baryons to be a simple perfect fluid 
in the curved spacetime (Equation ^). The stress- 
energy tensor describing such an fluid is given by 
( iMisner, Thorne fc Wheeler 1975| ) 



T^"' = pg^"' + {p + p)U^'U'' 



(5) 



where p and p are respectively the total proper energy 
density and pressure in the comoving frame. The is 
the four- velocity of the plasma fluid. The baryon- number 
and energy-momentum conservation laws are 



(6) 
(7) 



where ns is the baryon- number density. The radial com- 
ponent of Equation (|^) reduces to 



| + |(b + p)f/U) + ^|:(.^(p + p)c/'U) 
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The component of the energy-momentum conservation 
equation (|^) along a flow line is 



0(9) 



Equations (|^) and (g) give rise to the relativistic hydro- 
dynamic equations. 

We now turn to the analysis of e+e" pairs initially cre- 
ated in the Dyadosphere. Let be the proper densities 
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of electrons and positrons(e ). The rate equation for 
is 



{n,±U^),^ ^av[n,-{T)n,+ {T) 
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where a is the mean pair annihilation-creation cross- 
section, V is the thermal velocity of e*, and 7ie±(T) are 
the proper number-densities of e"^, given by appropriate 
Fermi integrals with zero chemical potential. The equilib- 
rium temperature T is determined by the thermalization 
processes occuring in the expanding plasma fluid with a 
total proper energy-density p, governed by the hydrody- 
namical equations (^j9[). We have 



Pt + Pe+ + Pe- + Pe- + PB, 



(11) 



where pj is the photon energy-density and Pe± is the 
proper energy-density of e^. In Equation (pj]), +Pb are 
baryon-matter contributions. We can also, analogously, 
evaluate the total pressure p. We define the total proper 
internal energy density e and the baryon mass density ps 
in the comoving frame, and have the equation of state (F 
is thermal index) 



e = p~PB, pB = nBmc, F 
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The calculation is initiated by depositing the total en- 
ergy (^ between the Reissner-Nordstrom radius r+ and 
the dyadosphere radius Vds- The calculation is continued 
as the plasma fluid expands, cools and the e'^e~ pairs re- 
combine, until it becomes optically thin: 



("-pairs + njardr = - 
R -J 
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where ctt is the Thomson cross-section, is the number- 
density of ionized electrons and integration is over the 
radial size of the expanding plasma fluid in the comoving 
frame. Here, we only present = 0, pB = case. 

We use a computer code (Wilson, Salmonson & Math- 
ews 1997, 1998) to evolve the spherically symmetric hy- 
drodynamic equations for the baryons, e"*'e~-pairs and 
photons deposited in the Dyadosphere. In addition, we 
use an analytical model to integrate the spherically sym- 
metric hydrodynamic equations with the following geome- 
tries of plasma fluid expansion: (i) spherical model: the 
radial component of four- velocity U{r) = U-^, where U is 
four- velocity at the surface {TZ) of the plasma, (ii) slab 1: 
U{r) = Ur = const., the constant width of expanding slab 
D = Ro in the coordinate frame of the plasma; (iii) slab 
2: the constant width of expanding slab R2 — Ri = Ro in 
the comoving frame of the plasma. 

We compute the relativistic Lorentz factor 7 of the ex- 
panding e+e^ pair and photon plasma. We compare these 
hydrodynamic calculations with simple models of the ex- 
pansion. In Figure (1) we see a comparison of the Lorentz 
factor of the expanding fluid as a function of radius for 
all of the models. We can see that the one-dimensional 
code (only a few significant points are pesented) matches 
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Fig. 1. Lorentz 7 as a function of radius. Three mod- 
els for the expansion pattern of the e^e~ pair plasma 
are compared with the results of the one dimensional hy- 
drodynamic code for a IOOOMq black hole with charge 
Q = O.lQjnax- The 1-D code has an expansion pattern 
that strongly resembles that of a shell with constant co- 
ordinate thickness. 



the expansion pattern of a shell of constant coordinate 
thickness (slab 1). 

We have shown that a relativistically expanding 
P.E.M. pulse can originate from the Dyadosphere of a 
EMBH. The P.E.M. pulse can produce gamma-ray bursts 
having the general characteristics of observed bursts. For 
example, the burst energy for a lOOOM© BH is 3 x 10^'* 
ergs with a spectral peak at 500 keV and a pulse dura- 
tion of 40 seconds (Ruffini, Salmonson, Wilson & Xue, 
1999). This oversimplified model is encouraging enough 
to demand further study of the Dyadosphere created by 
EMBHs. 
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